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I. INTRODUCTION 



Singlc-transvcrse-spin asymmetries (SSAs) have received much attention in recent years, both experimentally and 
theoretically. Large SSAs have been consistently observed in various experiments at different collision energies [H, H] • 
As a consequence of the parity and time-reversal invariance of the strong interaction, the SSA is directly connected to 
the transverse motion of partons inside a polarized hadron. Understanding the QCD dynamics behind the measured 
asymmetries should have the profound impact on our knowledge of strong interaction and hadron structure 4|. 

Two QCD mechanisms for generating SSAs have been proposed [HUH] and been applied [7l. M. Itj [ToL [Til. [l2Ul3l. Il4j 
extensively in phenomenological studies. Both of these mechanisms connect the SSA to the parton's transverse 
motion inside a transversely polarized hadron. One mechanism relics on the transverse momentum dependent (TMD) 
factorization for the observed polarized cross sections, and explicitly expresses the SSA in terms of "asymmetric" 
TMD parton distributions, known as the Sivers functions Q. The other follows the QCD collinear factorization 
approach for cross sections when all observed momentum scales are much larger than the non-perturbative hadronic 
scale 1/fm ~ Aqcd, and attributes the SSA to the twist-three transverse-spin dependent multi-parton correlation 
functions [1, 0j- Unlike the first mechanism, in which the SSA measures the TMD parton distribution and the spin- 
dependence of parton's transverse motion at a given momentum, the twist-three transverse-spin-dependent correlation 
functions reveal the net spin-dependence of parton's transverse motion when its transverse momentum is integrated. 
Naturally, the moment of th e sp in-dependent TMD parton distributions could be related to the twist-three multi- 
parton correlation functions [15| . Although two mechanisms each have their own kinematic domain of validity, they 
describe the same physics in the region where they overlap [l6l |. 

Most studies in both mechanisms have been concentrated on the SSAs generated by either the quark Sivers function 
[1, E(| [HI, EH EH or the twist-three quark-gluon correlation function @, H, E3, E3 , which is defined as 

Tf(x,x) = J ^^e^ p+ y^(p, s± \mi + b s±anfl F a + (y2)} ^(yi)\P,s±) , (i) 

with the gauge links suppressed. Possibilities of accessing the transverse motion of gluons, or, the gluon Sivers 
functions have also been investigated recently [H, EU . Likewise, the spin-dependence of the gluon's transverse motion 
in the QCD collinear factorization approach is represented by the twist-three tri-gluon correlation function, Tq(x, x), 
defined as 

T G (x, x) = J ^L^ e ^ + «r_l_(P ;Sx |^+ Q (0) [e s ^F+{y^)\ F a +(y^)\P, s ± ) . (2) 

Its contribution to the SSA was first studied by Ji in the context of direct photon production in hadronic collision [Tgj . 
Although direct-photon production provides a nice possibility to access Tg(x,x), the extraction of Tg{x,x) could be 
difficult due to the contribution from the quark-initiated subprocesses, and the limited knowledge on Tp{x, x). In this 
paper, we study the single-transverse spin asymmetry for open charm production in the semi-inclusive deep inelastic 
scattering (SIDIS) and argue that the SSA in SIDIS is a clean observable to extract the twist-three transverse-spin 
dependent tri-gluon correlation function, Tq{x,x). 

The _D-meson production at large enough transverse momentum P^± in SIDIS is dominated by the photon-gluon 
fusion subprocess, 7* + g — > c + 5, since the intrinsic charm contribution (20| is less relevant at large Ph± and the 
photon-charm subprocess, 7* + c — ► c + g, is suppressed by the small charm quark distribution at the collision energy 
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that we are interested in this paper. We calculate the SSA for D-meson production in lepton-proton SIDIS in terms of 
the QCD collinear factorization approach, and find that the asymmetry is directly proportional to the diagonal part 
of the twist-three tri-gluon correlation function in the polarized proton, Tq{x,x), because the photon-gluon fusion 
subprocess at this order does not have the so-called "hard-pole" contribution to the asymmetry Therefore, the 
measurement of SSA for D-meson production in SIDIS is a direct measurement of the tri-gluon correlation function, 
Tg{x, x). With a simple model for the tri-gluon correlation function, obtained under a similar assumption guiding the 
modeling of the quark-gluon correlation function, Tp(x, x), we find that the asymmetry for both COMPASS [2l[ and 
cRHIC [22| kinematics is sizable and could be measured experimentally. Recently, COMPASS experiment successfully 
measured the gluon polarization, AG, in a longitudinally polarized proton based on the photon-gluon fusion process 
by tagging charmed mesons [2l| . This certainly makes the measurement of SSA of open charm production in the 
same experimental setting and the extraction of the transverse-spin dependent tri-gluon correlation function, Tq(x, x), 
promising. 

We find that the SSA of D-meson production in SIDIS has a minimum at Zh ~ 0.5, and it increases as z% is moving 
away from this central value. This increase of the SSA away from Zh ~ 0.5 has the same physics origin as the observed 
increase of the magnitude of SSA in hadronic pion production as a function of increasing xf (or rapidity y), and is a 
prediction of the twist-three formalism of the QCD collinear factorization approach to the SSA. 

We also find that the twist-three gluonic contribution to the SSA has both similarity and difference from the 
twist-three fermionic contributions. Both gluonic and fermionic twist-three contributions to the SSA have the so- 
called "non-derivative" and "derivative" terms, which correspond to the terms that are proportional to the twist- 
three correlation functions and the derivative of the correlation functions, respectively. As noticed in Refs. [l7l. l23j. 
the fermionic "non-derivative" and "derivative" terms can be combined together and the dependence on the non- 
pcrturbative twist-three correlation function, Tp(x, x), is proportional to a simple combination, Tp(x, x) — xT' F (x, x). 
However, our explicit calculation shows that the "non-derivative" and "derivative" gluonic contribution can not be 
combined into the same simple form due to the difference in the partonic hard parts of these two terms. 

The same approach discussed in this paper can be applied for studying the SSA in open charm production in 
hadronic collisions, which is dominated by the gluon-gluon fusion if T@(x, x) is not too small comparing to the Tp(x, x) 
[HI US HH • With the extraction of the tri-gluon correlation function, Tg(x,x), and the knowledge of Tp(x,x), we 
enter a new era of exploring non-perturbative physics beyond the parton distribution functions (PDFs) which have 
been well-studied in the past thirty years. 

The rest of our paper is organized as follows. In Sec. |TT1 we present our calculation of the SSAs for open charm 
production in SIDIS. We first introduce the relevant kinematics of open charm production in SIDIS and present the 
formula for the unpolarizcd cross section. We then derive the twist-three formula for the SSA in QCD collinear 
factorization approach and express the asymmetry in terms of the tri-gluon transverse-spin dependent correlation 
function, Tq(x,x). We close this section by a discussion on the calculation of the color factor of the partonic hard 
part which depends on the color contraction of three gluon fields in the definition of the tri-gluon correlation function. 
In principle, there could be two gauge invariant tri-gluon correlation functions, Tq(x 7 x) and Tg(x,x) defined later, 
due to two independent ways to neutralize the color of the three gluon fields in the matrix elements. We point out that 
only one of them, Tq{x, x), could be related to the gluon Sivers function In Sec. IIII1 we estimate the production 
rate of open charm mesons in SIDIS for both COMPASS and eRHIC kinematics. We choose a simple ansatz for the 
tri-gluon correlation function, Tg(x,x), and present our predictions for the SSAs of open charm production at the 
existing COMPASS experiment and the planned eRHIC experiment. Finally, we conclude our paper in Sec. IIVI 



II. CALCULATION OF SINGLE-SPIN ASYMMETRY 



We start this section by specifying our notation and kinematics of SIDIS. We consider the scattering processes of 
an unpolarized lepton, e, on a polarized hadron, p, 

e(£) + p(P, a±) e(£') + h(P h ) + X, (3) 

where s± is the transverse spin vector defined below, h represents the observed D meson with momentum Ph and 
mass mh- For the collision energy that we are interested in this paper, we work in the approximation of one-photon 
exchange, and define the virtual photon momentum q = i — £' and its invariant mass Q 2 = —q 2 . We adopt the usual 
SIDIS variables: 

S ep = (P + £) , x B = — , y = —— = — , z h = — . 4 

2P ■ q P-i x B S ep P-q 
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It is also convenient to introduce the "transverse" component of the virtual photon momentum, q, as 

4t-Q p. p / p. Ph F h> W 
which is a space- like vector and orthogonal to both P and Ph . We dehnc 

1 q-P 



q± = -Qtitti = Q 



l 



x B P ■ Ph 



1L 

y-~ ' 



(6) 



To completely specify the kinematics, we will work in the so-called hadron frame [27] . where the virtual photon and 
the polarized proton are taken to have only one spatial component that is in the ^-direction: 

P" = P + »", tf = -XsP+n" + — (7) 

2x B P+ 

where the light-cone momentum component is defined as P^ = (P° ± P s )/\^2, and = (l + ,0~,0^), = 
(0 + , l _ ,0j_) are two light-like vectors with n ■ n — 1. The momentum of final-state D- meson can be written as 



j~.il XbP + 2 -a i z hQ a j~.il 

Ph = ^ mh±n + + P >^ (8) 



where rn\^_ = m,\ + Pfi ± with P^±_ = J P%j_. From Eq. ([6]) one can show that q± = \fqf_ = Phx/zh m this hadron 



h±- 

frame, independent of mass m/,. 

In this hadron frame, usually, one chooses the coordinate system such that the virtual photon has a vanishing 
energy component, corresponding to P + = Q/\/2xb 7 and Ph lies in the xz-plane (known as the hadron plane), as 
shown in Fig. [T] The lepton momenta, I and £' define the lepton plane and can be expressed in terms of variables ip 
and 4> as follows (27j . 

= ^ (cosh?/', sinh-i/j cos <f>, sinh ip sin <f>, — 1) , 
£ lpi = — (cosh?/', sm h ijj cos <j>, sinh ip sin (/>, +1) , (9) 



where </> is the azimuthal angle between the hadron and lepton plane, as indicated in Fig. [1] and 

coshV> = ^^-l = --l. (10) 

Q y 

We parametrize the transverse spin vector of the initial proton s± as 

s_l = (0, cos</> s ,sin0 s ,O), (11) 

where </> s is the azimuthal angle of sj_ measured from the hadron plane, as shown in Fig. [T] If one uses the lepton 
plane as the reference to define the azimuthal angle of s± as $5, and that of hadron plane as &h, one has the relation 
cj) s = $ s - $ h and 4> = -<5> h . 

The single transverse-spin asymmetry is defined as 

An = g ( s -L) - = dAa(s ± ) I da 

cr(s±) + <r(— sj_) dxsdydzhdPl^dcf) / dxsdydzhdP^dcp 

In the following subsections, we will first review the unpolarized cross section at leading order, and then derive the 
single-transverse polarized cross sections, Aa(s±). 

A. Unpolarized cross section 

The unpolarized differential SIDIS cross section may be calculated from the formula 

da Ka 2 em y 



dxBdydzhdPft ± d(/) Q 



ia L^(£,q)W^(P,q,P h ), (13) 
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where L^ v and are the leptonic and hadronic tensors, respectively. The leptonic tensor is given by 

g) = 2 (Vt + - ,9^Q 2 /2) ■ (14) 
The hadronic tensor has the following expression in QCD: 



W^(P, q, P h ) = ^E/ -0^e lq< (P\J»(O\XP h )(XP h \J u {Q)\P), 



(15) 



where is the quark electromagnetic current and X represents all other final-state hadrons other than the observed 
open charm meson h. 

The hadronic tensor can be further decomposed in terms of five parity and current conserving tensors V^ v [13] : 



(16) 



i=l 



where the Wi are structure functions which may be projected out from by W% — W p<r V?° ' , with the corresponding 
inverse tensors V,. Both V, and V, can be constructed from four orthonormal basis vectors: 



1 

Q 
l 



Zh 

' Z v X p T a , 



xbP" 1 



Y 1 ' 

Z» = (17) 

with normalization T 2 = 1 and X 2 = Y 2 = Z 2 = —1, which are reduced to those in [27| when rrih = 0. The tensor 
V5 does not contribute to the cross section when it is contracted with a symmetric L^, the other four tensors and 
their inverse are given as [27j : 



Vf v = X^X" + Y^Y\ = + Z»Z V , 

v »u = T ^ x u + T v XP ,^ yl xv = x ^ x u _ y^Y v , 

Vf = - (2T t "T u + X^X V + Y^Y V ) , V% v = T^T V 



y 3 



_1 (T»X V +T U X») 



l v = i (X^X U - Y^Y V ) 



(18) 



(19) 



The contraction of and Vf leads to various angular distributions. Let Ai = L llu V^ v / Q 2 , we have 



Ai = 1 + cosh 2 V, 



-4a 



cos</>sinh 2ip, A4 — cos 20 sinh ?/>. 



(20) 
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We can then write the cross section in Eq. JT3J) as 

da 



2 4 

dx B dydz h dP^d(j) Q 2 ^ 



AiWi. 



(21) 



At large Ph±_ ~ Q, the collinear factorization is expected to be valid, and Wi can be factorized into a convolution 
of the parton distribution function, the fragmentation function for the produced D meson, and a short-distance 
partonic hard part. The lowest-order (LO) contribution to the partonic hard part comes from the photon-gluon 
fusion subprocess 7* + g — ► Q(p c ) + Q(ps), which gives the leading order cross section as 



da 



dxsdydzhdP? 



dx 
x 



dz 



■G{x)D(z)6 



p2 



2 2 2 

— Q + z m r 



i\ 4 

- ^AW,, (22) 



where <7o — e c a em a sy/(8 7rz h'3 2 ): ^ = x b/x, z = Zh/z, and e c and m c are the fractional charge and mass of the charm 
quark, respectively. The P 2 ±/z 2 m the 5-function could be replaced by q\, and the 1/2 is the color factor. In Eq. |22|) 
G(x) is the unpolarized gluon distribution function with gluon momentum fraction x, and D(z) is the fragmentation 
function for the charm quark to become a D meson with z = P ■ Ph/P • Pc- We have suppressed the dependence on 
the factorization and renormalization scales for simplicity. We used Ph± ~ zp c ± inside the (^-function, which fixes the 
z integration. The lower limit of x integration x m i n is given by: 



x B 



XB 



z h (l-z h )Q- 
2m 2 1 



if z h + Jz 2 + -Jf > 1; 



p, 2 



p, 2 



iiz h + Jzi + ^<l 



(23) 



The short-distance parts Wi are calculated from the photon-gluon scattering and arc given by 



Wi 

m 
w 4 



u t 
- + - 

t u 

16x 2 " 



2sQ 2 



4i 2 s 



tv 



Ami 



4:xz—(u — i) 



3 -2 2 



t fi 



u t 
- + - 

t u 

tu 
1 

- + 

t 



where s,i,u are defined at the partonic level as 
1-x 



Q 2 - It Q 2 - 2u 



t 2 



= (xP + q) 2 = 



t={p c -q) 2 - m 2 c = -L-±Q 2 , 



2x^ 
Q 2 





i 






'1 


1" 


(!- 


- - - 




-8m* 




- - 




u 




t 


u 



(xP-p c ) 2 -m 2 = --C 
x 



(24) 



(25) 



which arc different from some definitions used in the literature. We found that this definition makes the expression 
of Wi for massive quark production sim pler . Taking m c = in Eqs. (j2"4")l and (|25[) . one recovers the results for the 
production of massless quark derived in |28l . I29I ] . 



B. Twist-three polarized cross section 



We now pro ceed to derive the single transverse-spin dependent cross section by applying the same method developed 
in Refs. [7lll7|. When both physically observed scales Q, Ph± Aqcd, the spin-dependent cross section for D-meson 
production is expected to be factorized in terms of twist-three transverse-spin dependent tri-gluon correlation function 



dAa(s±_) on 



2S e 



dzD(z) / dxidx27G(xi,X2) ie ps±nn lim — p-H(xi, X2, k±), 



(26) 
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where l/2S ep is the flux factor and € ps±nfi = e p<T ^ s^ a n^n u , 



P + dy 1 dy 2 ixiP + y - +i(x 2 -x 1 )P+ y - 
2n 



d aP {P,s ± \A a (0) [e s ^F+(y 2 )} A^)\P, 



(27) 



where d a p — —g a p + n a np + npn a . Tq{x\,X2) is related to the tri-gluon correlation function through Tq(x,x) = 
xTq(x, x). Since Tg{x\, X2) is real, we need an imaginary part of the hard-scattering function H(x\, x%, k±) to contract 
with i i P s ^ nn i n order to obtain a real Aa(s±). This imaginary part comes from the interference between a real part 
of scattering amplitude with a single gluon initial state and an imaginary part of the partonic scattering amplitude 
with an extra gluon, see Fig. [2] Technically, the imaginary part, or the phase, "i", arises when the virtual momentum 
integral of the extra gluon is evaluated by the residue of an unpinched pole from a propagator in the amplitude with 
an extra gluon. Such propagator is indicated by the one marked with a short bar in the diagrams in Fig. [3] 





Q V 








Q ' 











FIG. 2: A typical diagram that gives a non-vanishing contribution to the SSA. 

There are a total of eight partonic diagrams contributing to the twist-three polarized cross sections, Aer(s^). Four 
of them are shown in Fig. [3J and the other four are obtained by attaching the extra gluon in the same way on the 
right side of the final-state cut. When the extra gluon is attached to the left side of the final-state cut, as shown in 
Fig. [3 the phase from the propagator marked by the bar arises effectively as 



(p c - (%2 - x\)P - fcjj - ml + ie 2P • p c x\ - x 2 + v 1 ■ kj_ + ie 

■S(xi — X2+V1- k±), 



O(kl) 



2P-p, 



(28) 



to fix the virtual loop momentum fraction X\ — x 2 — V\ ■ k±_ with uj 1 = —2p p /2P ■ p c . On the other hand, the on-shell 



(xi -xi)p+h 



(a) 



(xi -xi)p+h 



s 



x 2 P+k L xiP 



(b) 




x 2 P+k L x,p 




(xi -xi)p+h 



(d) 



x 2 P+k± 



FIG. 3: Feynman diagrams that give the twist-three contribution to the spin-dependent cross section. The short bar indicates 
the propagator that produces the pole. The letters, A, B and C, represent the color of the initial-state gluons. 



condition associated with the unobserved anti-charm quark fixes the momentum fraction of the active initial-state 
gluon as 



8{pl-m 2 c ) = 5 ((x 2 P + k ± + q - pc) 2 - ml) 
= 2P-{ q -p c ) 6{X2 ~ X ' V2 - k ^ 



where terms at 0{k\) and higher are neglected and 

_ (q - Pc) 2 - ml 
X ~ 2P-(q-p c ) 



2p£ 
2P ■ (q - Pc 



(29) 



(30) 
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When the extra gluon is attached to the right hand side of the cut, the phase arises as 

1 11 



(p c + (%2 - x\)P + fcj.) - rnl - it 2P • p c x 2 - Xi - «i • fc_L - ie 

5(x 2 - x\ — v\- k±), 



+ 0(kl) 



2P-p, 

and the on-shell condition of the unobserved anti-charm quark gives 

1 



2P-(q- Pc ) 



5(x 1 - x), 



(31) 



(32) 



which has no fc^-dependence. 

Applying the so-called "master formula" in Ref. [l7j], we have from Eq. (|26[) the following general expression: 

lim ^7p" I dx Y [ dx 2 T G (xi,x 2 ) [H L (x 1 ,x 2 ,k 1 _)5(x 1 - x 2 + «i • k±)5(x 2 - x — v 2 ■ k±) 
k±^o dk'^ J J 

-HR(xi,x 2 ,k±)5(x 2 —xi — vi- kx)8(xi - x)} 

= [y2 - vi ) P H L (x, x, 0) — 
ax 



x lim 



d 



k ± ^0 dk P : 



[H L (x + (v 2 - Vi) ■ kj_,x + v 2 ■ k±, k±) - H R (x, x + vi- kj_,ki_)] , 



(33) 



where we have already used the facts that Hl{x, x, 0) = Hr(x, x, 0) and T G (x, x) = xT G (x, x). The fact that Eq. ([33)) 
depends only on the diagonal part of the tri-gluon correlation function, T G (x±, x 2 ), with x\ = x 2 = x is a consequence 
of that the photon-gluon fusion subprocess at this order has only the so-called "soft-pole" contribution to the SSA 
@; 13- Therefore, the measurement of the SSA in D- meson production in SIDIS is a direct measurement of the 
tri-gluon correlation function, T G {x,x). 

In terms of s, t, u defined in the previous subsection, we have 



< = — p£, 

u 



2x 
t 



■P'r 



(v 2 - fi) A 



2x 
i 



1 



P'c 



(34) 



Using Eqs. (|26[) . (|33[) . and adding the contributions from the eight diagrams together, we find the final expression for 
the fully differential singlc-transvcrsc-spin-depcndcnt cross section: 



dAa(sj_) 



dxBdydzhdPf 



= cr 



dx 



z 

V zi 



dz , N r / " 2 
D(z)5 



^- (1 '" )a ' £) Q 2 + £V 



1 + - 

u 



d ( T G (x,x)\ j (T G (x,x) 



dx V x 



\ x 



Ni 



, (35) 



where 1/4 is the color factor, T G (x, x) is the tri-gluon correlation function defined in Eq. ([2]), W% are given in Eq. I|24[). 
and the hard parts for the "non-derivative" term, Ni, are given by 



N 3 
N 4 



= 4 



2m 2 c ~ Q 2 6.t 2 



tii 









G- 


- - - 






u 





16x 2 









(! 


+ - - 


►')] 




u 





2Q a 

- — 



4z 2 g 2 _ 
tii 



2mt \ \ + - 



Q' 2 



t uj Q 



2zq\ 



tii 



+ s 



91 

iu 



1 1 

i u 



(36) 



Eq. (|35[) is our main result for the leading order twist-three T G (x, x) contribution to the fully differential polarized cross 
section, Aa(s±), of D-mcson production in SIDIS. The single transverse-spin asymmetry for the D-meson production 
in SIDIS is obtained by substituting Eqs. fl22J and (ESD into Eq. (fig) . 
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Similar to the twist-three contributions to the SSAs generated by the fermionic quark-gluon correlation function, 
Tp{x : x), the gluonic twist-three contribution to the SSA of D-meson production in Eq. ([55]) has both the "derivative" 
and "non-derivative" terms, a unique feature of twist-three contribution. It was found that the fermionic "non- 
derivative" and "derivative" terms can be combined into a simple form Tp(x,x) — xT F (x,x) (l7l. [23j. However, from 
Eqs. and (f3"u| . it is clear that Wi ^ Ni- That is, our explicit calculation shows that such a simple combination 
does not hold for contributions from tri-gluon correlation function Tq(x, x), and is not universal. 

We close this section by a discussion on the calculation of the color factor 1/4 in Eq. (|35j) . The color factor in 
Eq. (|35p depends on how colors of the three gluon fields in the matrix element of the tri-gluon correlation function 
in Eq. © are neutralized. If the color of the operator, F A (0) F c (y^) F B (y^), in Eq. @ with the Lorentz indices 
suppressed, is neutralized by {F c )ab = —if CAB with f CAB the fully antisymmetric structure constant of the color 
SU(3) group, we refer to this tri-gluon correlation function as Tq(x,x) as expressed in Eq. The corresponding 
color factor for the partonic part in Eq. (|3l>|) is calculated by contracting the color indices of the Feynman diagrams 
in Fig. [3] with N ^_ 1 s ) f ABC , which gives the color factor 1/4 in Eq. (f?5|) . On the other hand, if the color of the 

operator, F A (0) F c (y^) F B (lJi), is neutralized by (V c )ab = d CAB , which is symmetric under the interchange of any 
two indices, we have a different tri-gluon correlation function and we refer it as Tg{x, x), which has the same expression 
as Tq{x, x) except the difference in the contraction of the gluon color. The color factor for the corresponding partonic 
hard part is calculated by contracting the color indices of the same Feynman diagrams with ^ N i_ 1 ^ N i_ 4 } ) d ABC , which 
also gives the color factor 1/4. 

That is, there could be two tri-gluon correlation functions depending on how the colors of the three gluon fields 
are neutralized [3l| . We find that both correlation functions are gauge invariant after inserting necessary gauge link 
in the adjoint representation between the gluon field strengths in the matrix element in Eq. ©, and they contribute 
to the SSAs with the same partonic hard parts, and potentially, different color factors [26|]. Since the color factors 
are the same in our case, adding the potential contribution from Tq{x, x) is to replace the Tq{x, x) in Eq. (|35p by 
Tg(x,x) + T g (x,x). 

We noticed that all Feynman diagrams for producing a charm quark at this order in Fig. [3] have the same color 
structure, Tr[T A T c T B ] = (d ACB + i/ ACS )/4, which leads to the overall color factor 1/4 for the contributions from 
both Tg{x, x) and Tg(x, x). However, for the SSAs of producing a D meson, which is fragmented from an anticharm 
quark, both the partonic part and the antisymmetric part of the color structure change sign, while the symmetric 
part of the color structure is unchanged. Therefore, the SSAs for L>-meson production in the SIDIS at the leading 
order have the same functional form as that for the £>-meson produciton except the sum of the tri-gluon correlation 
functions, Tq{x,x) + Tq{x,x), is replaced by Tg{x,x) — Tg{x,x). That is, by comparing the SSAs for producing D 
and D mesons in SIDIS, we could gain valuable information on both tri-gluon correlation functions. However, the 
relation could be complicated in the £>-meson production in hadronic collisions due to the additional color How from 
the other colliding hadron [2(| . 

We also notice that the correlation function, Tq(x, x), with the color neutralized by (J- c )ab = —if CAB could 
potentially be related to the spin-dependent TMD gluon distribution, or the gluonic Sivers function [18|], since the 
middle gluon field strength of the operator, F A (0) F c \y^) F B (yi), could be related to the gauge link in the adjoint 
representation that is needed to define the TMD gluon distribution. Without knowing the size and sign of either 
tri-gluon correlation functions, we will treat them as one combined tri-gluon correlation function, labeled by Tq{x, x), 
in the following numerical estimates of the SSAs. 



III. PHENOMENOLOGY 



In this section we first evaluate the inclusive D- meson production rate at large Ph± in SIDIS. We then propose a 
simple model for tri-gluon correlation function Tq(x,x) and estimate the size of SSA for the D- meson production in 
SIDIS. 

The charm meson's transverse momentum, Phj_, is chosen to be along the x-direction in the hadron frame, and 
therefore, 

e PHS±nn = _p hxsin(j)s (37) 

The fully differential cross sections in Eqs. (f22|) and (f35|) can be decomposed in terms of the independent angular 
distributions as follows, 

dx B dydtdP^ = + * + 2*, 
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dAc 



sin 4> s (Acq + Aa\ cos <fi + A02 cos 2c, 



(38) 



Before evaluating the SSA, we first estimate the D-meson production rate in the unpolarized SIDIS by using our LO 
formula in Eq. (j2"2"]) . 

For the following numerical evaluations we use CTEQ6L parton distribution functions [32| . and charm-to- D frag- 
mentation functions from Ref. [33| ■ We choose the factorization scale to be equal to the renormalization scale and set 

1.3 GeV. In the following plots, we choose two sets of kinematic 



with charm quark mass m c 



variables. The first one is S ep = 300 GeV 2 , xb 
The other is S ep = 2500 GeV 2 , x B = 0.01 and Q 



0.01 and Q = 1 GeV, which is close to the COMPASS kinematics. 
= 4 GeV, which is more relevant to the planned eRHIC experiment 
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FIG. 4: The fully differential unpolarized cross section for D° production in SIDIS for COMPASS kinematics. The curves 
represent: og^solid), (j{ (dashed) , and ^(dotted) in Eq. (|38|) . 



In Fig. 21 we show individual coefficients of the angular distribution, <Tq , , and , of the fully differential 
unpolarized cross section for D° production in Eq. (|38[) as a function of both and Ph± for the kinematics relevant 
to COMPASS experiment. It is clear that the angular dependent pieces <§C <Jq , and might be too small to 

be significant. Without worrying about the detection efficiency, the D-meson production at Ph± ~ 1 GeV could 
be measurable. Likewise, Fig. O shows the fully differential unpolarized cross section for D° production for eRHIC 
kinematics. With larger Q and P/jj_, the production rate is smaller but may still have enough events with a high 
luminosity. 

In order to obtain the numerical estimate for the SSAs of D-meson production, we have to model the unknown, but 
universal, tri-gluon correlation function Tq(x, x). Similar to the ansatz for quark-gluon correlation function Tp(x, x), 
which was originally introduced in Q and found to be consistent with the latest experimental data [l7j . we model 
the tri-gluon correlation function Tq(x,x) as 



T G (x,x) = X g G(x) 



(39) 



with G(x) the normal unpolarized gluon distribution function. Because of its non-perturbativc nature, Tq{x, x) should 
be extracted from the experiments and the value and the sign of X g should be fixed by the data. For the following 
numerical estimate, we assume that X g has a positive sign and the same size as that for quark-gluon correlation 
function Tp(x,x) 0, and choose A g = 0.07 GeV. 

In order to present the SSA and its angular dependence on the </>, the angle between the hadron plane and the 
lepton plane, we define the ^-integrated single spin azimuthal asymmetries as 



(cos(n0)) 



1 lo^ cos ( n< ? 



dAa(s±) 



dx BdydzhdPr. d<p 



sm < 



Jo* d( t> 



da 



(40) 



dx Bdydzh,dP? ± c 



which gives 



(41) 
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In Fig.[5]we plot the SSAs as a function of Zh (left) and P/ l ± (right) for the COMPASS kinematics. The asymmetries, 
(1), (cos cf)), and (cos 20), defined in Eq. (jlTj) , are shown by the solid, dot-dashed, and dotted curves, respectively. 
For a comparison between the size of the "derivative" and the "non-derivative" terms, we also show, by the dashed 
curves, the contribution to the SSA, (1), from the derivative term only. It is clear that the derivative term dominates 
over the whole kinematic region. The asymmetries, (cost/)) and (cos 20), are too small to be observed experimentally. 
The SSA, (1), is of the order of 10%, and could be measurable at COMPASS experiment. 

Fig. [6] indicates that the SSA hits a minimum at ~ 0.5 and increases very fast when Zh becomes very large or 
very small. This is because the SSA, (1) ~ 1/(1 — x m i n ), due to the derivative of Tq(x,x) Q- From the definition 
of x m in in Eq. (|23[) . the Zh(l — Zh) has a maximum at Zh = 0.5. Therefore, x m m increases, equivalently, the SSA 
increases when Zh becomes either smaller or larger than 0.5. When Zh is much further away from the central value 0.5, 
the x m in becomes so large that the pcrturbatively calculated asymmetry could increase sharply, which could signal 
a breakdown of the twist-three approximation and a need of higher power corrections. Nevertheless, the increase of 
the SSA when z^ is moving away from the central value 0.5 has the same physics origin as the observed increase of 
the SSA as a function of increasing xf (or rapidity y) in the hadronic pion production [H, and it could be tested 
in the COMPASS experiment. 

Fig. [5] also indicates a monotonic increase of the SSA as a function of Phx- Although we expect the SSA to fall 
when Ph± increases, a natural behavior of the twist-three effect in QCD collinear factorization, the enhancement from 
the derivative of the Tq(x,x) at large x wins over the suppression from large Phi. due to the limited phase space at 
COMPASS kinematics. As we will see below, the decrease of the SSA as the increase of Ph± is clearly seen at the 
eRHIC kinematics. 

Similarly, we plot the SSAs for D° production for the eRHIC kinematics in Fig. [7] Due to the higher collision energy, 
the effective gluon momentum fraction x that dominates the SSAs is smaller, which leads to a smaller derivative of 
Tq(x, x) and a smaller SSAs. Similar feature has been seen in the SSA for hadronic pion production when we compare 
the data from the fixed-target experiments with that from RHIC experiments. The 5% SSA for D-meson production 
at eRHIC could be significant. 

The slightly different shape of the SSA as a function of Zh is purely a consequence of the difference in effective 
range of parton momentum fraction x. That is, the z/j-dcpcndcncc of the SSA provides a good measurement of the 
x-dependence of the correlation function, Tg(x, x). On the other hand, the slow falloff of the SSA as a function Phi. 
is natural due to the asymptotic \ g j Phi. behavior of the twist-3 contribution when Phi. increases. Of course, as 
discussed above, the 1/(1 — x m i n ) dependence of the twist-three formalism compensates some of the 1/P n ± falloff due 
to the phase space limit on parton momentum fraction x. 



IV. CONCLUSIONS 



In summary, we have studied the single transverse-spin asymmetry for _D-mcson production in SIDIS. In terms of 
QCD collinear factorization approach, we calculated both derivative and non-derivative contributions to the SSAs. 
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FIG. 6: Single-transverse-spin-asymmetries denned in Eq. (|41[l for D° production in SIDIS for COMPASS kinematics. The 
curves are: solid-(l), dashed-(l) with derivative-term only, dot-dashed- (cos <j>), and dotted- (cos 2(f)). 
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FIG. 7: Single-transverse-spin-asymmetries defined in Eq. (|41[) for D° production in SIDIS for eRHIC kinematics. The curves 
are: solid-(l), dashed-(l) with derivative-term only, dot-dashed- (cos cj>), and dotted- (cos 2(/>). 



At large enough transverse momentum, Ph±, the intrinsic charm contribution to the asymmetry might be neglected, 
and the SSA is directly proportional to the transverse-spin dependent tri-gluon correlation function, Tq{x,x) (or 
Tg(x,x) ± Tg(x,x) if we include both color structures), which has not been studied experimentally. We pointed 
out that by comparing the SSAs for producing D and D mesons in SIDIS, we could gain valuable information on 
both tri-gluon correlation functions. With a simple model for the Tg(x,x), we presented our estimates of the SSAs 
for the kinematics relevant for both COMPASS and future eRHIC experiments. From the inclusive production rate 
for the Z?-meson production and the estimated size of the SSAs, we argue that the SSAs of Z?-meson production in 
SIDIS could be a direct and clean probe of the unknown tri-gluon correlation function, Tq{x,x), which provides the 
important information on the spin-depenence of gluon's transverse motion inside a polarized hadron. 

However, we stress that the SSAs shown in all figures are directly proportional to the value and the sign of the 
X g and our model for the twist-three tri-gluon correlation function, Tg{x,x). A different x-dependence from that of 
G(x) could lead to a different derivative of Tq(x, x) and a different prediction for the SSA. The actual sign and size 
of the SSA, and the function, Tq[x,x), should be determined by the experimental measurements, just like the PDFs. 
However, our calculation does predict the short-distance dynamics and the kinematic dependence of the SSAs, such 
as the increase of the SSA when z% moves away from the central value 0.5. 

Finally, we emphasize that the QCD collinear factorization approach to the SSAs allows us to calculate the SSAs 
of open charm production or other particle production in hadronic collisions. With the experimental extraction of 
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the tri-gluon correlation function, Tg(x,x), as well as Tg(x,x), and the existing and new knowledge of Tp(x,x), we 
will be able to explore non-perturbative physics, in particular, the multi-parton quantum correlations, beyond what 
have learned from the parton distribution functions. 
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